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Abstract: Spatio-temporal assessment of the above ground biomass 
(AGB) is a cumbersome task due to the difficulties associated with the 
measurement of different tree parameters such as girth at breast height 
and height of trees. The present research was conducted in the campus of 
Birla Institute of Technology, Mesra, Ranchi, India, which is predomi¬ 
nantly covered by Sal (Shorea robusta C. F. Gaertn). Two methods of 
regression analysis was employed to determine the potential of remote 
sensing parameters with the AGB measured in the field such as linear 
regression analysis between the AGB and the individual bands, principal 
components (PCs) of the bands, vegetation indices (VI), and the PCs of 
the Vis respectively and multiple linear regression (MLR) analysis be¬ 
tween the AGB and all the variables in each category of data. From the 
linear regression analysis, it was found that only the NDVI exhibited 
regression coefficient value above 0.80 with the remaining parameters 
showing very low values. On the other hand, the MLR based analysis 
revealed significantly improved results as evidenced by the occurrence of 
very high correlation coefficient values of greater than 0.90 determined 
between the computed AGB from the MLR equations and field-estimated 
AGB thereby ascertaining their superiority in providing reliable estimates 
of AGB. The highest correlation coefficient of 0.99 is found with the 
MLR involving PCs of Vis. 

Keywords: above ground biomass; spectral response modeling; vegeta¬ 
tion indices; principal component analysis; linear and multiple regression 
analysis. 

Introduction 

Two types of remotely sensed data that have been used fre¬ 
quently by the researchers to extract information about biomass 
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and stand parameters comprise the optical data and non-optical 
data such as Radar and Lidar data. Roy et al. (1996) employed 
linear and multiple regression analysis techniques to determine 
biomass in Madhav National Park located in Madhya Pradesh, 
India from the spectral reflectance values of different bands of 
Landsat 5 TM data and the various vegetation indices and found 
that significantly high correlation exists between the biomass and 
middle infrared band and NDVI. On the other hand, among the 
various tasseled cap indices, greenness index and wetness index 
showed positive relationship with the biomass while the bright¬ 
ness index was found to be negatively correlated. Dong et al. 
(2003) also found positive relationship of NDVI with biomass. 
Lu (2006) described a host of techniques such as linear and non¬ 
linear regression, k-nearest neighbour, and neural network to 
determine the degree of relationship between biomass and spec¬ 
tral data. Le'vesque (1999) used fine resolution remotely sensed 
data such as Airborne digital camera and found that integration 
of spatial transforms of image fractions (texture, semivariogram 
range) and raw brightness values in high-resolution remotely 
sensed data have a good potentiality in modeling of forest struc¬ 
ture and health. Lefsky et al. (2001) evaluated the utility of sev¬ 
eral medium-resolution remotely sensed data for estimating stand 
structure attributes like age, basal area, biomass, and diameter at 
breast height (dbh). Their investigation revealed that multitempo¬ 
ral TM should be investigated as an alternative to either hyper- 
spectral or hyperspatial sensors, which are more expensive and 
more difficult to process than multitemporal Landsat TM. Only a 
few researches have been conducted regarding the biomass esti¬ 
mation using the coarse resolution data such as NOAA Ad¬ 
vanced Very High Resolution Radiometer (AVHRR), SPOT 
VEGETATION, and Moderate Resolution Imaging Spectro Ra¬ 
diometer (MODIS). However, coarse resolution data could be 
significantly useful for rapid estimation of biomass at continental 
and global scales (Fraser et al. 2002). 

The present study was carried out for the following objec¬ 
tives: To determine the coefficient of regression between the 
above ground biomass measured from the field with the spectral 
reflectance in individual bands of IRS LISS III data, individual 
principal components generated from the different spectral bands, 
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seven individual vegetation indices (Table 1) were considered in 
the present study, and the individual principal components gen¬ 
erated from seven vegetation indices respectively employing 
linear regression analysis; To develop multiple linear regression 
equations between the above ground biomass measured from the 
field with the spectral reflectance in different bands of IRS LISS 
III data, principal components generated from the different spec¬ 
tral bands, seven vegetation indices, and the principal compo¬ 
nents generated from seven vegetation indices respectively em¬ 
ploying multiple linear regression analysis; To determine the 
correlation coefficients between the above ground biomass com¬ 
puted through the multiple linear regression equations developed 
for four different types of remotely sensed parameters and the 
above ground biomass measured from the field in order to assess 
to what extent the various remote sensing parameters could pro¬ 
vide reliable estimates of the field-estimated above ground bio¬ 
mass. 


Materials and Methods 

Study area 

The study area selected for performing the present research cov¬ 
ers the campus of the Birla Institute of Technology Mesra, Ran¬ 
chi located within the Ranchi city, the capital of Jharkhand, India 
(Fig. 1) and is covered in the Survey of India toposheet No. 73 
E/7. The study area represents natural and protected open type 
forest with Sal as the dominant species. The study area is sur¬ 
rounded by Subarnarekha and Jumar River. There are three main 
seasons that prevail in the region: summer season (15 March - 
15 June) when the climate is extreme with maximum tempera¬ 
ture of 38°C and minimum 21°C, rainy season (June - October) 
with the annual rainfall being about 1530 mm and the winter 
season (November - February) with temperature ranging between 
2°C and 22°C. The vegetation remains green during the rainy 
and winter season while the summer season is characterized by 
both leaf shedding and leaf emerging phases. The soil type with¬ 
in the study area is largely lateritic (Kumar 1999). 

Materials used 

To perform the present investigation, LISS III data of IRS-P6 
satellite (Path/Row- 105/55) acquired on 20th October, 2008 was 
used. The LISS III data comprises four spectral bands viz. Green 
(0.52-0.59 pm), Red (0.62-0.68 pm), NIR (0.77-0.86 pm) and 
MIR (1.55-1.70pm) each having 23.5m spatial resolution and 8- 
bit radiometric resolution. Erdas Imagine 9.1® was employed for 
performing all digital image processing tasks such as georefer- 
encing of the raw image, subsetting of the plots, extraction of 
reflectance values, generation of vegetation indices and principal 
component images. The ARCMAP 9.2® was used to generate 
canopy density and species type maps. The geographic coordi¬ 
nates of the different sampling plots selected within the study 
area were determined using Leica Handheld Global Positioning 
System. The linear regression and multiple regression analyses 
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between the biomass and the other governing parameters were 
performed using the SPSS 17.0 software. MS-Excel software 
was used for generation of plots between the biomass and vege¬ 
tation parameters and remote sensing variables such as reflec¬ 
tance values in different bands, vegetation indices and PCs of 
bands and PCs of vegetation indices. The original satellite image 
was georeferenced with the GCPs extracted from the Survey of 
India toposheet No. 73E/7 on 1:50 000 scale. 
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Fig. 1 Study area showing the locations of sampling plots 

Methods used 

The task of biomass estimation of the Sal trees through remote 
sensing involved two broad steps as mentioned below. 

Methods used for field data collection and above ground bio¬ 
mass estimation 

Collection of various field data pertaining to the Sal trees and 
subsequently, their above ground biomass estimation involved 
the following steps. 

Keeping in view the extent of the study area and the variation 
in the canopy density of the Sal trees within the study area, ten 
representative sampling plots were selected. At the time of se¬ 
lecting the plots it was kept in mind that they should be large 
enough to ensure that there was no influence by edges, road, 
and/or pipelines. 

The size of each sampling plot was taken as 1 hectare, which 
is approximately equal to 16 pixels of the LISS III data. Exact 
geographic coordinates of these plots in the field were deter¬ 
mined using GPS (Leica Handheld Global Positioning System). 

Total number of trees existing within each sampling plot was 
counted. 

The gbh (girth at breast height) of the individual trees occur¬ 
ring within each plot was measured using measuring tape and 
was later converted to dbh (diameter at breast height) by dividing 
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the gbh values by n (0.318). 

An important parameter for the volume estimation of trees 
through volumetric equation is the height of tree. The measure¬ 
ment of the tree height was carried out through indirect approach. 
First, using a hypsometer the viewing angle between the observer 
and position at the first forcking of the tree was measured and 
then the distance between the observer and the base of the tree 
was measured using a measuring tape. The height of the tree was 
determined using the following formula. 

Tan 0 = height of the tree / distance from tree (1) 

where 0 denotes the viewing angle. 

The volume of each tree was determined using the following 
formula, 

F = 0.118 + 0.257 D 2 *H (2) 

where, V is the volume of Sal trees, D is the diameter at breast 
height, and H is the height at the first forcking. 

The total above ground biomass of a plot is calculated as the 
sum of the total bole biomass of all the trees occurring within the 
plot and the litter biomass contained in the plot. The bole bio¬ 
mass represents the dry weight of the standing tree whereas the 
litter biomass comprises the dry weight of the leaves lying on the 
ground. The bole biomass of each tree was determined by multi¬ 
plying the volume of the tree with its specific gravity. For Sal 


trees the specific gravity has been determined as 0.726 (Shrestha 
2003). The litter biomass within each plot was determined using 
the following procedure. The litter lying within a few small sam¬ 
pling plots of size (lmxlm) selected at different locations within 
the main plot were collected from the ground and then their dry 
weight was measured in the laboratory which was then averaged 
to obtain the litter biomass of a lxlm plot. The total litter bio¬ 
mass of the main plot was determined by multiplying the average 
litter biomass corresponding to lxlm plot with the area of the 
main plot. Then the total bole biomass and total litter biomass 
estimated for each main plot was summed up to determine the 
total above ground biomass within each main plot. 

Methods employed for estimation of above ground biomass using 
the parameters derived from remotely sensed data 
The satellite data was georeferenced with the help of the Survey 
of India toposheet No. 73 E/7 on 1:50 000 scale with sub-pixel 
accuracy. Then each sampling plot selected in the field was cor¬ 
rectly delineated on the satellite image using the coordinates 
measured in the field through GPS. 

Different vegetation indices such as NDVI, SAVI, RVI, DVI, 
IPVI, GNDVI and GVI (Table 1) were generated through model 
maker tool in ERDAS IMAGINE 9.1 software. This was fol¬ 
lowed by the generation of principal components of the bands 
and stacked vegetation indices respectively. 


Table 1: Equations of different vegetation indices used in this study 


Vegetation Indices 


Explanation of symbols 

References 

NDVI 

(NIR - R) 

~ (NIR + R) 


NDVI=Normalized Difference Vegetation Index, NIR = near 

infrared, R = red, G = green, and L = Soil-adjustment factor 

Rouse et al. (1974) 

RVI = 

NIR 


RVI (Ratio Vegetation Index) 

Richardson and Wiegand (1977) 

R 


SAVI = 

(NIR - R) 

+ (l + L ) 

SAVI (Soil-adjusted Vegetation Index) 

Huete (1988) 

~ (NIR + R + L) 

IPVI = 

NIR 


IPVI (Infrared Percentage Vegetation Index) 

Crippen (1990) 

(NIR + R) 


GNDVI 

_ (NIR - G) 

~~ (NIR + G) 


GNDVI (Green Normalized Vegetation Index) 

Jensen (2006) 

GVI = 

(NIR - R) 

-+ 128 

GVI (Green Vegetation Index) 

Jensen (2006) 

(NIR + R) x 127 

DVI = 

(NIR - R ) 


DVI (Difference Vegetation Index) 

Lillesand et al. (2007) 


Then the average value of the reflectance in different bands, 
different vegetation indices, individual PCs generated from the 
bands and stacked vegetation indices for each sampling plot were 
computed. 

Both linear and multiple linear regression analyses were per¬ 
formed between the average AGB value and different remote 
sensing parameters as mentioned in the preceding step for the 
individual plots. Regression plots were also generated between 
the AGB and individual remote sensing parameters. Multiple 
linear regression analyses were performed employing the SPSS 
17.0 software. 

Multiple linear regression equations were developed between 


the field-estimated AGB for each of the ten different plots se¬ 
lected in the entire study area and the corresponding remote sens¬ 
ing based parameters viz. the four spectral bands of LISS III data, 
four principal components generated from the spectral bands, 
seven vegetation indices considered in the present study, and the 
seven principal components determined from the seven vegeta¬ 
tion indices respectively. Then AGB was computed for each 
sampling plot by substituting the corresponding reflectance val¬ 
ues of four spectral bands, PC values of the spectral bands, seven 
vegetation indices and the PC values of the seven vegetation 
indices in the respective multiple linear regression equations. 
Correlation analysis was performed between the field-estimated 
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AGB and the AGB estimated through the multiple linear regres¬ 
sion analysis for different sampling plots to determine to what 
extent the regression based AGB tallies with the field estimated 
AGB and consequently, to determine to what extent the remotely 
sensed parameters could provide correct estimates of the AGB 
measured in the field. 

Finally, comparative analysis was performed among the re¬ 
sults obtained from the various regression analyses in order to 
determine the particular spectral band, vegetation index, and 
PC(s) of bands and vegetation indices that exhibited the highest 
regression coefficient with the field-estimated AGB. 

Results 

The results determined from the present research have been pre¬ 
sented in two sections; section one comprises the discussion 
about the linear regression analysis performed between AGB and 
individual remote sensing parameters and section two deals with 
the multiple linear regression analysis. 

Linear Regression Analysis 

As mentioned earlier, linear regression analysis was performed 
between the field-estimated AGB determined from each of the 
ten sampling plots and the corresponding remotely sensed pa¬ 
rameters such as the average reflectance values recorded in dif¬ 
ferent spectral bands, individual vegetation indices, individual 
PCs generated from the spectral bands and the vegetation indices 
respectively (Table 2). Comparative analysis of the regression 
coefficient values (R 2 ) given in table indicates the occurrence of 
significantly higher R 2 values with the vegetation indices and 
their PCs as compared to that exhibited by the different spectral 
bands and their PCs. Among the four spectral bands of LISS III 
data used in the present study, Red band showed the highest 
regression coefficient with the field-estimated AGB with R 2 val¬ 
ue of 0.65 while the remaining three bands exhibited nearly same 
R 2 values (Fig. 2). Similarly, the highest value of regression coef¬ 
ficient (0.66) is exhibited by PC2 while the remaining three PCs 
show nearly similar R 2 values ranging between 0.56 and 0.59 
(Fig. 3). On the other hand, the different vegetation indices ex¬ 
hibit significantly much higher regression coefficients with the 
AGB as compared to the spectral bands and the PCs of the spec¬ 
tral bands. Among the different vegetation indices used in the 
study, NDVI shows the highest regression coefficient (0.81) 
followed in decreasing order by SAVI (R 2 = 0.78), IPVI (R 2 = 
0.75), GVI (R 2 = 0.73), RVI (R 2 = 0.71), DVI (R 2 = 0.70) and 
GNDVI (R 2 = 0.66) (Fig. 4). The PCs of the seven vegetation 
indices also exhibit considerably high regression coefficient with 
the AGB which is found to be lower than that exhibited by the 
vegetation indices but significantly higher than the regression 
coefficient values exhibited by the spectral bands and the PCs of 
the spectral bands respectively. The first PC of the vegetation 
indices exhibits the highest regression coefficient with a value of 
0.77 followed in decreasing order by PC2 (R 2 = 0.72), PC3 (R 2 = 
0.67), PC4 (R 2 = 0.64), PC6 (R 2 = 0.55), PC5 (R 2 = 0.54) and 
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PC7 (R 2 = 0.52) respectively (Fig. 5). 


Table 2: Regression coefficient values of different 
parameters with field-estimated AGB 

remote sensing 

Remote Sensing Parameters 

Individual Parameters 

R 2 


Green 

0.60 


Red 

0.65 

Spectral Bands 


NIR 

0.62 


MIR 

0.59 


PC 1 

0.57 


PC 2 

0.66 

PCs of Spectral Bands 


PC 3 

0.59 


PC 4 

0.56 


NDVI 

0.81 


SAVI 

0.78 


RVI 

0.71 

Vegetation Indices 

DVI 

0.70 


IP VI 

0.75 


GNDVI 

0.66 


GVI 

0.73 


PC 1 

0.77 


PC 2 

0.72 


PC 3 

0.67 

PCs of Vegetation Indices 

PC 4 

0.64 


PC 5 

0.54 


PC 6 

0.55 


PC 7 

0.52 



Fig. 2 Relationships between AGB and spectral bands 



Fig. 3 Relationships between AGB and PCs of spectral bands 
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Fig. 4 Relationships between AGB and vegetation indices 



Fig. 5 Relationships between AGB and PCs of vegetation indices 

From the foregoing analyses it can be inferred that the field- 
estimated AGB is best regressed with the NDVI values while the 
SAVI also exhibited significantly high regression coefficient 
with the AGB. The other remote sensing parameters such as the 
spectral bands, PCs of the spectral bands, remaining five vegeta¬ 
tion indices and PCs of the seven vegetation indices failed to 
provide considerable regression coefficients with the AGB. 
Therefore, NDVI and SAVI can be considered as reliable estima¬ 
tors of the AGB (Majumdar 2010). 

Multiple Linear Regression Analysis 

Multiple linear regressions were developed between the AGB 
estimated for each of the ten different sampling plots and the 
corresponding average DN values of each of the four remote 
sensing based parameters considered in the present study viz. the 
spectral bands, PCs of the spectral bands, seven vegetation indi¬ 
ces, and the PCs of the seven vegetation indices and are shown in 
equations (3), (4), (5) and (6) respectively. 

AGBgands = 0.006* Green + 0.051* Red + 0.043* NIR + 0.020* 
MIR-4.707 (3) 

AGBpcs of Bands = 0.018* PC 1 + 0.019* PC 2 - 0.007* PC 3 + 
0.009* PC 4-6.218 (4) 

AGBvegetation Indices = 0.024* NDVI - 0.027* SAVI + 0.009* RVI 
+ 0.004* DVI + 0.024* IPVI + 0.001* GNDVI + 0.004* GVI - 
7.663 (5) 


AGBpcs of Vegetation Indices “ 0.015* PC 1 + 0.002* PC 2 + 0.007* 
PC 3 + 0.011 * PC 4 + 0.006* PC 5 + 0.000* PC 6 - 0.001 * PC 

7 - 5.295 (6) 

Subsequently, the AGB of each sampling plot was computed 
by substituting the corresponding average DN values of the dif¬ 
ferent remote sensing parameters viz. the spectral bands, PCs of 
the spectral bands, seven vegetation indices, and PCs of the 
vegetation indices in the respective multiple linear regression 
equations. This is referred to as the computed AGB of the differ¬ 
ent sampling plots. The primary objective of the present research 
is to assess the potential of different remote sensing derived pa¬ 
rameters for providing the correct estimates of the AGB meas¬ 
ured from the field. In order to perform this task, the correlation 
coefficient values between the computed AGB derived from 
different multiple linear regressions and the AGB estimated from 
the field were determined (Figs. 6a, 6b, 6c, 6d). Comparative 
analysis of these figures revealed the occurrence of the highest 
correlation coefficient of 0.99 between the computed AGB based 
on the principal components of the vegetation indices and the 
field-estimated AGB (Fig. 6a), followed in decreasing order by 
the vegetation indices (0.96), spectral bands (0.92) and principal 
components of the spectral bands (0.91) (Figs. 6b, 6c, 6d respec¬ 
tively). These findings prove that the multiple linear regression 
derived on the basis of the PCs of the seven vegetation indices 
selected for the present investigation is capable of providing 
estimates of the field-measured AGB up to 100 percent accuracy 
while the multiple linear regressions derived based on the re¬ 
maining three remote sensing parameters could also be consid¬ 
ered as potential reliable estimators of the above ground biomass. 
The occurrence of the very high correlation coefficients of the 
computed AGB with the field-estimated AGB based on the mul¬ 
tiple linear regressions using all the four types of remote sensing 
data could be attributed to the combined contribution of the spec¬ 
tral information from the multivariable data and also, the preva¬ 
lence of monotonous vegetation (Sal) with nearly homogeneous 
canopy density thereby resulting in similar spectral characteris¬ 
tics throughout the study area. 


a: PCs of Vegetation Indices 


Fig. 6 Correlation between computed AGB and field-estimated AGB 
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Discussions 

The above ground biomass estimation from the remote sensing 
data performed in the present study involved two types of regres¬ 
sion analysis such as the linear regression analysis and multiple 
linear regression analysis. Based on the results of the linear re¬ 
gression analysis performed between the above ground biomass 
and different remote sensing parameters such as the individual 
spectral bands of IRS LISS III data, individual PCs generated 
from the spectral bands, seven vegetation indices, and the indi¬ 
vidual principal components generated from the seven vegetation 
indices, it was found that it is only the NDVI that exhibited a 
regression coefficient value of greater than 0.80 with the above 
ground biomass. In addition, the PCI of the vegetation indices 
showed a regression coefficient of 0.72. The other remote sens¬ 
ing parameters exhibited significantly very low regression coef¬ 
ficient values with the AGB. 

On the other hand, the multiple linear regression analysis per¬ 
formed between the field-estimated above ground biomass and 
various types of remote sensing parameters such as the spectral 
bands, PCs of the spectral bands, vegetation indices and PCs of 
the vegetation indices exhibited significantly higher values of 
regression coefficients. The correlation coefficient values deter¬ 
mined between the computed above ground biomass derived 
from the multiple linear regression equations and the field- 
estimated above ground biomass were found to be significantly 
very high; for example, 0.99 for the principal components of the 
vegetation indices, 0.96 for the vegetation indices, 0.92 for the 
spectral bands and 0.91 for the principal components of the spec¬ 
tral bands. The occurrence of significantly high values of correla¬ 
tion coefficients can be attributed to the combined contribution 
of the spectral information of the vegetation canopy from the 
multiple variables. Two salient findings that emerged from the 
present investigation are as follows. First, the multiple linear 
regression based analysis possesses better potential of providing 
correct estimates of the above ground biomass as compared to 
the linear regression based analysis and second, the principal 
components of the vegetation indices are better correlated with 
the above ground biomass. 
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